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Universal computation of a quantum system consisting of superpositions of well-separated coherent states 
of multiple harmonic oscillators can be achieved by three families of adiabatic holonomic gates. The first gate 
consists of moving a coherent state around a closed path in phase space, resulting in a relative Berry phase 
between that state and the other states. The second gate consists of “colliding” two coherent states of the same 
oscillator, resulting in coherent population transfer between them. The third gate is an effective controlled-phase 
gate on coherent states of two different oscillators. Such gates should be realizable via reservoir engineering of 
systems which support tunable nonlinearities, such as trapped ions and circuit QED. 
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Reservoir engineering schemes continue to reveal promis¬ 
ing new directions in the search for potentially robust 
and readily realizable quantum memory platforms. Such 
schemes are often described by Lindbladians [1] possessing 
decoherence-free subspaces (DFSs) [2] or (more generally) 
noiseless subsystems (NSs) [3] - multidimensional spaces im¬ 
mune to the nonunitary effects of the Lindbladian and, po¬ 
tentially, to other error channels [4, 5]. On the other hand, 
holonomic quantum computation (HQC) [6] is a promising 
framework for achieving noise-resistant quantum computa¬ 
tion [7]. In HQC, states undergo adiabatic closed-loop parallel 
transport in parameter space, acquiring Berry phases or matri¬ 
ces (also called non-Abelian holonomies or Wilson loops [8]) 
which can be combined to achieve universal computation. 

It is natural to consider combining the above two concepts. 
After the initial proposals [9, 10], the idea of HQC on a DFS 
gained traction in Refs. [11,51] and numerous investigations 
into HQC on DFSs [12] and NSs [13-15] followed. How¬ 
ever, previous proposals perform HQC on DFS states con¬ 
structed out of a finite-dimensional basis of atomic or spin 
states. There has been little investigation [16] of HQC on 
DFSs consisting of nontrivial oscillator states (e.g. coher¬ 
ent states [17, 18]). While this is likely due to a historically 
higher degree of control of spin systems, recent experimen¬ 
tal progress in control of microwave cavities [19-21], trapped 
ions [22], and Rydberg atoms [23] suggests that oscillator- 
type systems are also within reach. In this Letter, we propose 
an oscillator HQC-on-DFS scheme using cat-codes. 

Cat-codes are quantum memories for coherent-state quan¬ 
tum information processing [24] storing information in super¬ 
positions of well-separated coherent states which are evenly 
distributed around the origin of phase space. Cat-code quan¬ 
tum information can be protected from cavity dephasing via 
passive quantum error correction [25] using Lindbladian- 



Figure 1. In the d = 2 cat-code, quantum information is encoded 
in the coherent states |tr 0 (0)) = \a) and |oq(0)) = \-a). (a) Wigner 
function sketch of the state before (top) and after (bottom) a loop 
gate acting on |—or), depicting the path of \-a) during the gate (blue) 
and a shift in the fringes between | ±a). (b) Phase space diagram for 
the loop gate; X = \{a + cf) and P - -\{a - o'). The parameter 
aft) is varied along a closed path (blue) of area A, after which the 
state |-or) gains a phase 6 - 2A relative to \a). (c) Effective Bloch 
sphere of the | ±a) qubit depicting the rotation caused by the d = 2 
loop gate. Black arrow depicts initial state while red arrow is the 
state after application of the gate. The dotted blue arrow does not 
represent the path traveled since the states leave the logical space 
| ±a) during the gate, (d-f) Analogous descriptions of the collision 
gate, which consists of reducing a to 0, driving back to a exp (if), 
and rotating back to a. 


based reservoir engineering [5]. In addition, such information 
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can be actively protected from photon loss events [5, 26-28]. 
While there exist plenty of methods to create and manipu¬ 
late the necessary states [5, 26, 29-31] and while the gates 
can also be implemented using Hamiltonians, we consider 
reservoir engineering due to its protective features. Cat-codes 
differ from the well-known Gottesman-Kitaev-Preskill (GKP) 
encoding scheme [32] in both state structure and protection. 
GKP codes consist of superpositions of highly squeezed states 
and focus on protecting against small shifts in oscillator po¬ 
sition and momentum. In contrast, cat codes protect against 
damping and dephasing errors, the dominant loss mechanisms 
for most cavity systems. While realistic GKP realization 
schemes remain scarce [33], cat-codes benefit from greater 
near-term experimental feasibility [21]. 

For simplicity, let us introduce our framework using a sin¬ 
gle oscillator (or mode). Consider the Lindbladian 

1 d ~ l 

p = FpF^ --{F^F,p} with F= (1) 

2 y=0 

[a,o'] = 1, h = a f a, k e R, dimensionless a y e C, and p a 
density matrix. The d = 1 case [F = s[k{cl - oro)] reduces to 
the well-known driven damped harmonic oscillator ([34], Sec. 
9.1) whose unique steady state is the coherent state |oro) (with 
a\ao) = ao|<ro». Variants of the d = 2 case are manifest in 
driven 2-photon absorption ([35], Sec. 13.2.2), the degener¬ 
ate parametric oscillator ([36], Eq. 12.10), and a laser-driven 
trapped ion ([37], Fig. 2d; see also [38]). A motivation for this 
work has been the recent realization of the F = yfi<(a 2 - afy 
process in circuit QED [20], following an earlier proposal to 
realize F = s[k( a d - a d ) with d = 2,4 [5]. For arbitrary d 
and certain a v , a qudit steady state space is spanned by the 
d well-separated coherent states \a v ) that are annihilated by 
F. The main conclusion of this work is that universal con¬ 
trol of this qudit can be done via two simple gate families, 
loop gates and collision gates , that rely on adiabatic varia¬ 
tion of the parameters a v (t). Universal computation on mul¬ 
tiple modes can then be achieved with the help of an entan¬ 
gling two-oscillator infinity gate. We first sketch the d = 2 
case and extend to arbitrary d with \a v ) arranged in a circle in 
phase space. The straightforward generalization to arbitrary 
arrangements of | a v ) is presented in [39]. We then discuss gate 
errors and integration with cat-code error correction schemes 
[5, 28], concluding with a discussion of experimental imple¬ 
mentation. 

Single-qubit gates. —Let d = 2 and let ao , oq depend on 
time in Eq. (1), so the steady-state space holds a qubit worth 
of information. The positions of the qubit’s two states | a v (t)) 
in phase space are each controlled by a tunable parameter. We 
let oro(0) = -oq(0) = a (with a real unless stated otherwise). 
This system’s steady states | ±a) are the starting point of pa¬ 
rameter space evolution for this section and the qubit defined 
by them (for large enough a) is shown in Fig. la. 

The loop gate involves an adiabatic variation of a fit) 
through a closed path in phase space (see Fig. lb). The state 
|a fit)) will follow the path and, as long as the path is well 


separated from \ao(t)) = | a), will pick up a phase 6 = 2A, 
with A being the area enclosed by the path [40]. It should be 
clear that initializing the qubit in \-a) will produce only an 
irrelevant overall phase upon application of the gate (similar 
to the d = 1 case). However, once the qubit is initialized in a 
superposition of the two coherent states with coefficients c±, 
the gate will impart a relative phase: 

c+\a) + C- | -a) —> c+\a) + C-e l ° \—a) . (2) 


Hence, if we pick | a) to be the x-axis of the | ±a) qubit 
Bloch sphere, this gate can be thought of as a rotation around 
that axis (depicted blue in Fig. lc). Similarly, adiabatically 
traversing a closed and isolated path with the other state pa¬ 
rameter |ofo(0) induce a phase on \a). 

We now introduce the remaining Bloch sphere components 
of the cat-code qubit. For a = 0, the d = 2 case retains its 
qubit steady-state space, which now consists of Fock states 
\p), p - 0,1 (since F - y[Ka 2 annihilates both). One may 
have noticed that both states \±a) go to |0) in the a —> 0 limit 
and do not reproduce the a = 0 steady state basis. This issue 
is resolved by introducing the cat state basis [41] 
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(k*> + i-f i-qt» 
(3) 


with normalization - y^[l + (-)^exp(-2 a 2 )]. As a —> 
0> | p a ) ~ I p) while for a —> oo, the cat states (exponentially) 
quickly become “macroscopic” superpositions of | ±a). This 
problem thus has only two distinct parameter regimes', one 
in which coherent states come together (a 1) and one in 
which they are well-separated (a » 1, or more practically 
a > 2 for d = 2). Eq. (3) shows that (for large enough a) cat 
states and coherent states become conjugate z- and x-bases 
respectively, forming a qubit. We note that p = 0,1 labels the 
respective ±1 eigenspace of the parity operator exp(/7ra); this 
photon parity is preserved during the collision gate. 

We utilize the a 1 regime to perform rotations around 
the Bloch sphere z-axis (Fig. If), which effectively induce a 
collision and population transfer between | a) and | -a). The 
procedure hinges on the following observation: applying a 
bosonic rotation R$ = exp (i(pn) to well-separated coherent 
or cat state superpositions does not induce state-dependent 
phases while applying R$ to Fock state superpositions does. 
Only one tunable parameter ao(t) = -a fit) is necessary here, 
so F = yfl<[d 2 - afit) 2 ] with |or 0 (0)| = a. The collision gate 
consists of reducing a to 0, driving back to a exp(/0), and ro¬ 
tating back to a (Fig. le). The full gate is thus represented by 
R^S^S^, with S(p [42] denoting the nonunitary driving from 0 
to a exp (i(f>). Since 


4Vo = R l(K* s oW = SoR l S l ’ (4) 


the collision gate is equivalent to reducing a , applying R J on 
the steady-state basis | p), and driving back to a. The net result 
is thus a relative phase between the states \p a ): 


co|O a > + ci|l a )—* c 0 |0 a ) + c\e ,0 |l a ). 


(5) 
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In the coherent state basis, this translates to a coherent popu¬ 
lation transfer between \±a). 


with yu = 0,1, • • • ,d - 1. The corresponding cat-state basis 
generalizes Eq. (3) to 




Figure 2. Sketch of the adiabatic paths of the components (a) | -Sa) 
and (b) \-S~ l a) during the infinity gate. 


Two-qubit gates .—Now let’s add a second mode b and in¬ 
troduce the entangling infinity gate for the 2-photon case. We 
now use two jump operators for Eq. (1), 

Ff = (a - a) (a + Sa) and Fu - (< ab-a 2 )(ab + Sa 2 ). (6) 


We keep a > 0 constant and vary S(t) in a figure-eight or “oo” 
pattern (Fig. 2), starting and ending with 6=1. For S = 1, the 
four DFS basis elements {|±or)}<8){|±ar)} are annihilated by both 
Fj and F n . For 6 ± 1 and for sufficiently large a, the basis 
elements become | a, a), \a, -Sa), \-6a, a), and | -6a, -S~ l a). 
Notice that the S~ l makes sure that F n \ - Sa, -S~ x a) = 0. 
This S~ l allows the fourth state to gain a Berry phase distinct 
from the other three states. Since Berry phases of different 
modes add, we analyze the a/b -mode contributions individu¬ 
ally. For any state which contains the | -6a) component (in ei¬ 
ther mode), the Berry phase gained for each of the two circles 
is proportional to their areas. Since the oppositely oriented 
circles have the same area (Fig. 2a), these phases will cancel. 
The Berry phase of the fourth state, which contains the com¬ 
ponent | -S~ l a), will be proportional to the total area enclosed 
by the path made by S~ l . Inversion maps circles to circles, 
but the two inverted circles will now have different areas (Fig. 
2b). Summing the Berry phases if/i gained upon traversal of 
the two circles i e { 1,2} yields an effective phase gate: 

| -a, -a) + |rest) —> e l( ^ ] \-a, -a) + |rest), (7) 

where |rest) is the unaffected superposition of the remaining 
components {| a, a), \a, -a), \-a, a)}. 

Single-qudit gates .—We now outline the system and its 
single-mode gates for arbitrary d. Here we let a y (0) = ae v 
with real non-negative a, e v = exp(/^v), and v = 0,1, • • • ,d— 
1 (see Fig. 3a for d = 3). This choice of initial qudit con¬ 
figuration makes Eq. (1) invariant under the discrete rotation 
exp (Ojn) and is a bosonic analogue of a particle on a discrete 
ring [43]. Therefore, hmodd is a good quantum number and 
we can distinguish eigenspaces of exp(z^fi) by projections 
[44] 


oo | d -1 

= F ' [dn+ ^ dn +a*i = -j F exp 

n =0 a v=0 


dn 

i—{n~n)v 

a 


( 8 ) 


n^la) ( \d) a -> 0 (9a) 

lMa> = ~ { Td ^ V \ae v ) a -» oo . (9b) 

Since overlap between coherent states decays exponentially 
with a, the quantum Fourier transform between coherent 
states | ae v ) and cat states \/i a ) in Eq. (9b) is valid in the well- 
separated regime, i.e., when 2a sin ^ » 1 (satisfied when 
\(a\aei)\ 2 <^c 1). It should be clear that the more coherent 
states there are (larger d), the more one has to drive to resolve 
them (larger a). Also note the proper convergence to Fock 
states \jd) as a —^ 0 in Eq. (9a). 




Figure 3. (a) Three-fold symmetric configuration of steady states |or v ) 
of Eq. (1) with d = 3 and depiction of a loop gate (6) acting on \a 2 ) 
and a collision gate (0) between all states, (b) Arbitrary configuration 
of steady states for d -1 , depicting |ar 0 ) undergoing a loop gate and 
|ori), \a 2 ) undergoing a displaced collision gate (see Supplementary 
Material). 


Both gates generalize straightforwardly (see Fig. 3a for d = 
3). The loop gate consists of adiabatic evolution of a specific 
a y (t) around a closed path isolated from all other ay (0). There 
are d such possible evolutions, each imparting a phase on its 
respective | ae v ). The collision gate is performed as follows: 
starting with the \ae y ) configuration for large enough a, tune 
a to zero (or close to zero), pump back to a different phase 
<rexp(/0), and rotate back to the initial configuration. Each 
| Ha) will g a i n a phase proportional to its mean photon number, 
which behaves in the two parameter regimes as follows: 

. ... . / ld + 0(a 2d ) o' —> 0 (10a) 

Wa\n\Ha) = ) , , 2 

{ a 2 + 0(a 2 e ) a -> oo , (10b) 

where c = 1 - cos 2i: ld. Since a rotation imparts only a yu- 
independent (i.e. overall) phase in the well-separated regime 
[Eq. (10b)], the only //-dependent (i.e. nontrivial) contribution 
of the symmetric collision gate path to the Berry matrix is 
at a = 0. This gate therefore effectively applies the Berry 
matrix exp (—iQx) to the qudit, where x = F\Fa)(Fa\ I s 
the discrete position operator of a particle on a discrete ring 
[43]. More generally, one does not have to tune a all the way 
to zero to achieve similar gates - e.g. being in the regime with 
2 a sin ^ « 1 is sufficient. The two-mode infinity gate can 
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likewise be extended to the d-photon case and it is a simple 
exercise to prove universality [39]. 

Gate errors. — In the Lindbladian-dominated adiabatic 
limit [39], the role of the excitation gap is played by the dis¬ 
sipation gap - the eigenvalue of the Lindblad operator whose 
real part is closest (but not equal) to zero. Since our Lind- 
bladians are infinite-dimensional, it is possible for the dissi¬ 
pation gap to approach zero for sufficiently large |ar v | (i.e., in 
the limit of an infinite-dimensional space). For the symmet¬ 
ric d-photon case however, this is not the case and the gap 
actually increases with a (verified numerically for d < 10). 
Having numerically verified the infinity-gate, we also see that 
the gap increases with a in the two-mode system (6). The gap 
can also be seen to increase by analyzing the excitation gap of 
the Hamiltonian F*F (see Ref. [15], Sec. VIII.C). 

Here we discuss the scaling of leading-order nonadiabatic 
errors, focusing on the single-mode gates for d = 2, 3. 

Non-adiabatic corrections in Lindbladians are in general 
nonunitary, so their effect is manifest in the impurity of the 
final state (assuming a pure initial state). Extensive numerical 
simulations [45] show that the impurity can be fit to 

e= l-Tr{p(D 2 }cx -J— (11) 

KTaP 

as a, T — > oo, where p(T) is the state after completion of the 
gate, a = |or v (0)| is the initial distance of all \a v ) from the 
origin, p > 0 is gate-dependent, and k is the overall rate of 
Eq. (1). One can see that e « 0(T ~ l ), as expected for a nearly 
adiabatic process. Additionally, we report that p « 1.8 for 
d = 2 and p « 3.9 for d = 3 loop gates, respectively. For the 
d = 2,3 collision gates, we observe that p « 0. 

Photon loss errors .—We have determined that the above 
gates can be made compatible with a (photon number) parity- 
based scheme protecting against photon loss [5, 28]. In such 
a scheme, one encodes quantum information in a logical 
space spanned by even parity states (e.g. | a v ) + \-a v ) with 
v — 0,1, • ■ ■ , d — 1, generalizing Sec. II.D.3 of [25]). Pho¬ 
ton loss events can be detected by quantum non-demolition 
measurements of the parity operator (-1)". In the case of 
fixed-parity cat-codes, errors due to photon loss events can 
be corrected immediately [28] or tracked in parallel with the 
computation [5]. By doubling the size d of the DFS of Lind- 
bladian (1) to accomodate both even and odd parity logical 
spaces, we have determined a set of holonomic gates which 
are parity conserving and are universal on each parity sub¬ 
space [39] . This allows for parity detection to be performed 
before/after HQC. 

Implementation & conclusion. —We show how to achieve 
universal computation of an arbitrary configuration of multi- 
mode well-separated coherent states | a v ) by adiabatic closed- 
loop variation of a v (t). We construct Lindbladians which ad¬ 
mit a decoherence-free subspace consisting of such states and 
whose jump operators consist of lowering operators of the 
modes. One can obtain the desired jump operators by non- 
linearly coupling the multi-mode system to auxiliary modes 
(c, d, • • •), which act as effective thermal reservoirs for the ac¬ 


tive modes. For the case of one active mode a , if one assumes 
a coupling of the form a$ + H.c. and no thermal fluctuations 
in c, one will obtain (in the Born-Markov approximation) a 
Lindbladian with jump operator a. Therefore, a generalization 
of the coupling to F$ + H.c. will result in the desired single¬ 
mode Lindbladian (1) with jump operator F. Since F are poly¬ 
nomials in the lowering operators of the active modes, quartic 
and higher mode interactions need to be engineered. Such 
terms can be obtained by driving an atom in a harmonic trap 
with multiple lasers [37] or by coupling between a Josephson 
junction and a microwave cavity [5, 20]. We thus describe ar¬ 
guably the first approach to achieve holonomic quantum con¬ 
trol of realistic continuous variable systems. 
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Summary 

In Sec. A, we provide a qudit extension of the two-qubit infinity gate described in the main text. In Sec. B, we prove 
universal quantum computation on the subspace spanned by coherent states lying on a circle in phase space using loop, collision, 
and infinity gates. In Sec. C, we relax the restriction of the coherent states lying on a circle and extend our gates to arbitrary 
configurations. In Sec. D, we make contact with the adiabatic theorem and Berry connections, revealing the Berry matrix for our 
gates. In Sec. E, we show how to implement gates of fixed photon number parity, allowing for photon loss errors to be corrected 
before/after the gates. 


A. Two-qudit gates 

Here we extend the two-mode oo gate to arbitrary d. The initial adiabatic path is the same, but an extra single mode loop gate 
is required afterwards to make the infinity gate a true controlled phase gate. 

Similar to the single-mode gate extensions stated in the main text, we let a v ( 0) = ae v be the initial coherent state values for 
each mode [with real non-negative a , e v = exp (/^v), and v = 0,1, ■ • • ,d - 1]. The two-mode jump operators from Eq. (6) of 
the main text now generalize to 


d-2 d -3 

Fj - (d - Sad- 1 ) J”[(^ _ &v) and F n = (db - aad-i)(ab - Saad-i) ]~[{db - aa y ). 

v=0 v=0 


( 1 ) 


When S = 1 , the d 2 DFS basis elements {\a v }} <S> {\a v Q} are annihilated by both Fj and F u . For 5 ± 1 and for sufficiently large 
a, the basis elements become {\a v , Xa v Q] and {\Sad-u Yoy)}, where X follows from Tab. I. Upon adiabatic evolution of S in an 
oo-shaped path, the Berry phase obtained from any \Sa v ) vanishes due to each of the two loops contributing an equal area with 
opposite sign. The Berry phase of the non-trivial terms, which contain the component \S~ l a v ), per loop i e { 1,2} is now 


xfi - 2 a 2 



( 2 ) 


which can be calculated using the corresponding Berry connections Tg s = /(d _1 Qjd|6| ?arg £|c) _1 (r) (see Appx. D). This calcula¬ 
tion is equivalent to the geometrical reasoning provided in the main text. At the end of the gate, the phase gained from adiabatic 
evolution of all d 2 DFS states is listed in Tab. II. We can see that all DFS states of the form \aj-i, a v ) gain the phase if/\ + if/ 2 , 
with exception of the last state \aj~i , ®d-\)- If we now apply a loop gate on the |oy-i) state of the first oscillator with a 
phase, the last state becomes the only one to gain a phase: 


\a d - U a d -i) + |rest> —> e \a d -i,a d - X ) + |rest>, 

where |rest) is the unaffected superposition of the remaining components {| a v , ay)}. 


( 3 ) 


B. Proof of universal computation 

To prove universal computation of the three effective unitary gate families, we analyze their generators and show that various 
commutators of those generators provide a basis for the Lie algebra su (d) [1]. The loop gates are generated by d coherent 
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Table I. Values of X for each state of mode b for all d 2 steady states of a Lindbladian with jump operators Fj and F n from Eq. (1). 
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Table II. Berry phase gained for each of the d 2 DFS states upon completion of the oo-gate. 


state projections n y = \ae v ){ae v \ and can be commuted with^ to generate su (d) as follows. Superpositions of ft v provide d - 1 
linearly independent traceless diagonal elements. The linearly independent off-diagonal coherent state transitions are provided 
by \ d(d - 1) elements 


8w' — [^V? — 


| ae v ){ae v | 
e~ i2 f(y~y') _ i 


+ H.c. 


(4) 


along with the \d(d - 1) elements ^\tc v - ft v ,g vV ] (both sets for v ^ V). Together this makes d 2 - 1 elements, the dimension of 
su(d). Adding in the qudit infinity gates from above then becomes sufficient for universal computation [2]. 




Figure 3 from main text, (a) Three-fold symmetric configuration of steady states \a v ) of Eq. (1) with d - 3 and depiction of a loop gate (6) 
acting on \a 2 ) and a collision gate (0) between all states, (b) Arbitrary configuration of steady states for d - 7, depicting |or 0 ) undergoing a 
loop gate and \a\), \a 2 ) (with “center of mass” y and relative coordinate a) undergoing a displaced collision gate. 


C. Generalization to arbitrary a v 

Here we discuss generalized versions of the single-mode symmetric 2-photon and d-photon gates from the main text, breaking 
symmetry and allowing for arbitrary values of |ck v ) (as long as all states remain well-separated). Fig. 3 from the main text is 
reproduced for convenience. 
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Generalized loop gates can be performed via simultaneous variation of multiple a v (t). The procedure is exactly the same as 
that described in the main text for the d = 2 loop gate. More interestingly, collision gates away from the origin between any 
two (or more) coherent states can also be performed. For example, the symmetric d - 2 collision gate, demonstrated with the 
collision at the origin in Fig. Id (see main text), can be performed anywhere in phase space. If one picks two states |oq(0)), 1^(0)) 
(see Fig. 2b) and assumes that all other states are well-separated from them, a two-state collision gate is then a displaced version 
of the original collision gate from Fig. Id. To better understand this, first let a\ = y + a and (%2 = y - a with complex relative 
coordinate a and “center of mass” y. Then, note that the displaced cat states become (up to a global phase) 

D y \n a ) '“^ 2 2=(|ai> + k* 2 », 

V2 

where /i e {0, 1} and D y is the displacement operator (with D y |0> = |y». The states D y \ji a ) provide a basis (valid for any aq, ar 2 ) 
for the steady states of the d = 2 case F = V ^(a - a i )(a - a 2 ) from Eq. (1) of the main text, which reduces to F = a 2 - a 2 ) 
when y = 0. The collision gate thus translates into varying the amplitude and phase of a in the path shown in Fig. 2b while 
keeping y constant, producing the same net result as the original d = 2 collision gate and inducing transitions between |oq(0)) 
and |or 2 (0)>. Naturally, such phase and collision gates allow for universal control for |or v (0)) at arbitrary locations in phase space 
(Fig. 2b), granted that |or v (0)) are well-separated from each other: |or v (0) - ay (0)| > 4 for all v, V - 0, 1, • • • ,d- 1. 


D. Berry connections 

The adiabatic theorem states that traversal of a closed loop in parameter space in the adiabatic limit (parallel transport) results 
in the application of a Berry matrix (more generally, a holonomy) on the initial state space, generalizing the Berry phase for 
non-degenerate systems. There are two adiabatic limits associated with Lindbladians: one governed by the steady states of the 
entire Lindbladian ([3], Sec. 6; see also [4-10]) and one governed by eigenstates of a Hamiltonian part [11]. We consider a 
DFS in the former limit. The Berry matrix imparted on a d-dimensional steady state subspace when time evolution is governed 
by a Lindbladian may not be equal to that imparted on the same subspace when time evolution is governed by a Hamiltonian. 
However, Ref. [6] showed that the Berry matrix of the Lindbladian-controlled DFS case [5] does indeed reduce to the ordinary 
SU(d) Berry matrix for unitary systems [12]. Here we sketch another proof of this result, based on a more general derivation of 
the Berry matrix for generic steady state subspaces from Ref. [10]. We note that non-adiabatic HQC schemes on DFSs [13] and 
NSs [14] have also been proposed and that HQC (as well as general manipulations) on a DFS can alternatively be understood in 
terms of a generalized quantum Zeno effect [6, 15-18]. 

The Lindbladian Berry matrix, in general nonunitary, can be calculated via an ordered path integral of the Lindbladian Berry 
connections 


pn' 'i (T(T ' 




(5) 


where is the conserved quantity [19] corresponding to \ji a )<ji' a \ and d\ = The determine the initial Bloch vector 
of the cat qudit and are known in closed form only for d = 2 [15]. Interestingly, do not participate in parallel transport 
for either gate and Eq. (5) reduces to its closed system counterpart - a superposition of ordinary Berry connections between cat 
states. This reduction holds because do not “cross-talk,” i.e., P a J^P^ = 0 with P a = YufZo \Pa)(Pa\ the projection on the 
steady state subspace and P^ = 1 - P a . This condition can be shown to imply that adiabatic parallel transport on the DFS is free 
from the nonunitary effects of the Lindbladian. We have shown the condition numerically for the example here; the exact proof 
is shown in Ref. [10]. In other words, one can verify that 


<2p 

‘yu// ;crcr' 


— X , , 4^ -I- X A A 

— U l±'(Teller ^ U H<T^n'tT' 


( 6 ) 


with = i(fi a \d^cr a ) the ordinary Berry connection [12]. 

The Berry connections A A reduce exactly to the gate generators^, ft v from the main text. For example, the symmetric collision 
gate in Figs, lc and 2a arises from changes in the magnitude and phase of a (here complex). Therefore, A e {| a\, arg a = (p) and 
a simple calculation reveals A^ f = 0 and 




-<V</*al %«> 


H 1^1 * 0 

-<VM 2 N -» 00 • 


(7a) 

(7b) 


One can see that A* oc% as \a\ —* 0, confirming that^ indeed generates the effective operation induced by the collision gate. 
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E. Integration with photon loss error correction scheme 


Here we discuss the parity conserving holonomic single-qubit gates; the infinity gate extends analogously. To add HQC 
capability to qudits protected from photon loss, we need to double the size of the DFS in order to accommodate both even and 
odd parity logical spaces. This is done by letting a - a v —> a 2 - in Eq. (1) from the main text. The original {|or v >} basis (for 
sufficiently large and well-separated a y ) becomes 


I«V, p) ~ (|av> + (-) p I-av» , (8) 

where p € {0,1} (mod 2) indexes the logical space and 2d is the dimension of the new DFS. In order to extend HQC to these states, 
we need to show that both gate families act identically on both logical spaces p. The loop gate generalizes straightforwardly. 
Varying ay (t) in a closed loop far away from all other a v produces the same geometric phase for both |±ay), so | ay,/?) — > 
e l 11 ° |ay, p) with 6 independent of p. To generalize the collision gate [for which we now set a v = a exp(z|v)], we once again need 
a dual basis to properly take the a —> 0 limit. This cat state basis {\p a , p) = |(2 p + p ) a )} is obtained by letting p —> 2p + p and 
d —> 2 d in Eq. (9) from the main text and plugging in Eq. (8). In the two regimes of interest, 


( \2p + p) ct —> 0 

I Zv^o e_ifv(2/i+p) l ae l >,p) a —> oo . 


(9a) 

(9b) 


The phase gained during the a = 0 rotation part of the collision gate is then cp(p a , p\h\p a , p) “~° cp(2p + p). Therefore, the 
collision gate also acts the same way on both logical spaces, up to an overall phase of exp (-icpp). 
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